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We give an estimate of the gravitational field of force exerted on a test particle by the far galaxies, 
in the frame of the weak field approximation. In virtue of Hubble's law, the action of the far matter 
turns out to be non negligible, and even the dominant one. A nonvanishing contribution is obtained 
only if the discrete and fractal nature of the matter distribution is taken into account. The force 
per unit mass acting on a test particle is found to be of the order of 0.2 cHq, where c is the speed 
of light and Ho the present value of Hubble's constant. 
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The main idea underlying the present paper comes 
from a critical reading of the first work in which the ex- 
istence of dark matter was proposed, namely the work of 
Zwicky |1] in which the virial theorem was applied to the 
Coma cluster. The point concerns the role of the forces 
due to the external matter. Let us recall that in the first 
use of the virial theorem, namely that of Clausius for 
a gas enclosed in a vessel, the main contribution to the 
virial was coming from the external forces (the pressure 
of the walls), because the contribution of the internal 
forces was either vanishing (for a perfect gas) or negligi- 
ble. In the case of Coma, the contribution of the internal 
forces was still found by Zwicky to be negligible, but the 
contribution of the external galaxies was neither taken 
into account, nor even mentioned at all, perhaps because 
one could not conceive how it may act as a pressure. The 
missing contribution to the virial was thus attributed to 
some hypothetical interior dark matter. In the present 
letter we show how it is conceivable that the external 
matter (actually, the far one) may produce a pressure, 
and we also give a rough estimate of its contribution to 
the virial, which appears to be in rather good agreement 
with the observations. 

The first point we make is that, as the forces on a test 
particle depend not only on the positions of the galaxies, 
but also on their velocities and accelerations, due to Hub- 
ble's law the dominating contribution comes from the far 
matter. This comes about as follows. From the point of 
view of general relativity, in the weak-field approxima- 



tion the problem of estimating the force (per unit mass) 
on a test particle amounts to writing down the equa- 
tions for the geodesic motion when the metric tensor g^^^, 
is a solution of the Einstein equation with the external 
galaxies as sources. Writing the metric tensor as a per- 
turbation of the Lorentzian background rj^^, , namely, as 
S/jf = VtJ.i' + hfii^j the perturbation h^j^^, turns out to be 
a solution of the wave equation, so that its components 
are analogous to the familiar retarded potentials of elec- 
trodynamics (although relevant differences exist between 
the two cases, as particularly emphasized by Zeldovich 
and Novikov [3| ). In fact one finds 
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where G is the gravitational constant, while Mj, 
and 7j , j — 1, . . . , N , are the mass, the position vector 
and the Lorentz factor of the j-th source galaxy, the dot 
denoting derivative with respect to proper time along the 
worldline of the source. 

As hfii, depends on each source not only through its 
position, but also through its velocity, the latter has to 
be assigned in order that the model be defined. To this 
end we use Hubble's law as a phcnomenological prescrip- 
tion. Making reference to a local chart with Lorentzian 
coordinates having as origin the center of mass of the 
considered localized system (the Coma cluster), for the 
velocities = qO) of the external galaxies we thus re- 
quire 



•Electronic address: 'carati@mat.ummi.itj 

t Electronic address: sergio.cacciatori@uni nsubria.it| 

■f Electronic address: jgalgani@mat.unimi.itl 



7, ^Hoqj 



3 = 1, 



(2) 



(the dot denoting now derivative with respect to the 
background Lorentzian time). Here, Hq is the Hubble 
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constant which, in our extremely simplified model, we 
take fixed at its present value. 

Notice that the Hubble assumption ([2]) has an essen- 
tial impact on the size of the gravitational field of force. 
Indeed such a force contains a term (decreasing as 
proportional to the velocity of the source, and a term 
(decreasing as 1/r) proportional to the acceleration of 
the source. Thus, estimating the acceleration through 
Hubble's law, the latter term actually doesn't depend on 
distance at all, so that the far matter is found to give 
the dominant contribution to the gravitational field of 
force. This situation is reminiscent of the way in which 
Mach's principle was dealt with in [3| (see page 102). The 
main difference being that in such a case, lacking Hub- 
ble's law, the velocities of the sources were neglected. 
Consequently, only the Newtonian, fast decaying, poten- 
tial was considered, so that only the near matter, and 
not the far one, appeared to play a role. 

So we address our attention to the dominating term, 
proportional to the acceleration of the source. Such a 
term, which we denote by f , has the form 

where the masses of the galaxies were all put equal to a 
common value M, and the Lorentz factors jj were put 
equal to 1, for the reasons to be illustrated later. No- 
tice the extremely simple nature of this force per unit 
mass (or acceleration). Apart from a multiplicative fac- 
tor, such a force is just the sum of all the unit vectors 
pointing to each of the external galaxies. Actually, our 
attention was addressed to the component of such a force 
f along a given direction. Such a component will be sim- 
ply denoted by /, and the corresponding component of 
u by u. 

Having determined the quantity of interest (/ or u), we 
come now to our second main point, namely the problem 
of how to describe the distribution of the external galax- 
ies. It is immediately seen that u exactly vanishes (at 
any point) if the external matter is described as a contin- 
uous medium with a spherically symmetric density. This 
should be expected, in agreement with BirkhofF's theo- 
rem. So we take a different point of view, analogous to 
the probabilistic one introduced by Chandrasekhar and 
von Neumann (see the review 0]) in the problem of es- 
timating the vector sum of the Newtonian forces exerted 
on a star by the surrounding ones. In such an approach, 
the external sources are conceived as point particles, and 
the place of the matter density of the continuum case is 
now taken by a probability density for the position of a 
galaxy. From such a point of view, the previous result 
(the vanishing of u for a spherically symmetric matter 
density) , now reads as the vanishing of the mean value of 
u for a spherically symmetric probability density of the 
position of a galaxy. 

We thus come to an estimate of the variance of the 
force / (or of u). It will be seen that the result depends 
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FIG. 1: Histogram of the random variable u, which is pro- 
portional to the force per unit mass due to A'' external galax- 
ies. The histogram was computed with 10, 000 samples of 
= 512, 000 galaxies, by counting the fraction of times the 
value u belongs to a given interval of width \/]V. 



on the further assumptions one introduces concerning the 
spatial distribution of the external galaxies. Assume first 
that the positions q-,- of the N galaxies are independent 
random variables, uniformly distributed with respect to 
the Lebesgue measure. Then the sum defining u is found 
to grow as -y/iV. This indeed is just a consequence of 
the central limit theorem, because in such a case u is 
the sum of N independent identically distributed ran- 
dom variables which turn out to have zero mean and a 
finite variance. By the way, such a result is the analogue 
of that obtained by Chandrasekhar and von Neumann 
for the case of Newtonian forces, the only difference be- 
ing that in their case the variance is infinite (due to the 
divergence of Newton's force at zero distance, a property 
which plays no role in our case). For what concerns the 
corresponding estimate of the virial, one easily sees that 
with the present assumption the estimate is by far too 
small to account for the observations, just because the 
considered sum behaves as Vn rather than as N (see 
later). 

So we modify the previous assumption, and consider 
the case in which the probability density is fractal Q 
(see also @). This means first of all that the positions 
of the galaxies are no more independently distributed, 
so that u is no more constrained to grow as \/N, and 
can instead have a faster growth, as required by the ob- 
servations. However, the analytical computation of the 
probability distribution of the field of force now becomes 
a quite nontrivial task, with respect to the much simpler 
case considered by Chandrasekhar and von Neumann. 
So we are forced, at least provisionally, to investigate the 
problem by numerical methods. 

We proceeded as follows. In order to estimate (at the 
origin of the coordinates) the sum defining u, the posi- 
tions of the N galaxies are extracted (with the method 
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FIG. 2: The variance of u versus the number A'^ of galaxies 
in log-log scale. The dashed line is the curve cr^ — 0.2 7V^. 



described in in such a way that the mass distribu- 
tion has a fractal dimension, precisely the fractal dimen- 
sion 2. The corresponding histogram, obtained through 
10,000 samples of configurations of the galaxies, is shown 
in Fig. [2 for N = 512,000. 

We then study the dependence of u on the number N of 
external galaxies, which was made to vary in the range 
1000 < < 512, 000, with the density kept constant. 
This means that the positions of the N points were taken 
to lie inside a cutoff sphere whose volume was made to 
increase as N. For the values of N investigated, the cor- 
responding radius turns out to be so small with respect 
to the present horizon, that the Lorentz factors 7 could 
altogether be put equal to 1 (as was previously assumed) , 
and more in general the special relativistic character of 
our model was actually justified. 

The mean of u turns out to practically vanish for all 
TV, while its variance cr^ is found to grow as N"^ (actually, 
as 0.2 A^^), rather than as as occurs in the uniform 
case. This is shown in Fig. [51 We thus conclude that the 
standard deviation u/ of the component of the force per 
unit mass along a direction is proportional to N , being 
given by 



MN = 



MN . 



(4) 



(where — c/Hq is the present horizon). 

We now take such a result, which was obtained for 
extremely small values of and extrapolate it up to 
the present horizon Rq — c/Hq^ i,e., we insert in formula 
(m the actual value of N, so that the quantity MN can 
be identified with the total visible mass of the Universe. 

Concerning the total visible mass MN of the Universe, 
one can write 
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shown (see later) that a quite natural consistency condi- 
tion of our model leads to the estimate 
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where po = (3Hq)/ (SttG), with JIq ~ 0.05, is the ac- 
tual density. Inserting this in ^ one gets a/ ~ 0.2 cHq. 

On the other hand, if a random variable / has zero 
mean and a finite variance with great probability its 
modulus will take on values very near to its standard 
deviation <7f. In such a sense we may say to have found 
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which perhaps constitutes the main result of the present 
work. Namely, in our oversimplified model within the 
fractal hypothesis, the force per unit mass, i.e., the ac- 
celeration, exerted by the far matter on a test particle, 
is found to have a value of the order of cHq, which is the 
one that is needed in most cases in which the presence of 
a dark matter is advocated. Notice that the assumption 
of a uniform, rather than fractal, distribution of mat- 
ter would lead instead to |/| ~ cHq/Vn, i.e., essentially 
to / ~ 0. Namely, without the fractal hypothesis the 
Zwicky procedure of neglecting at all the gravitational 
contribution of the external matter, would be justified. 

We can now apply our estimate to the case of the virial 
theorem for a cluster of galaxies. To this end, we have 
first of all to assume a property that might in principle 
be checked, namely, that the forces at sufficiently sepa- 
rated points are uncorrelated. It is then conceivable that 
locally, in some regions, such a random field of force may 
form patterns of a central type, attractive towards a cen- 
ter. This is obviously equivalent to admit that locally, in 
such special regions, the external far matter produces a 
pressure. 

Let us recall that, according to the virial theorem, for 
a confined system of n particles one has 
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here, — {1/n) J^i'^i variance of the veloc- 

ity distribution of the galaxies of the cluster, whereas 
V = ■ is called the virial of the forces (per unit 

mass) , Xj denoting the position vector of the i-th internal 
particle with respect to the center of mass of the cluster, 
and fj the force per unit mass acting on it, while overline 
denotes time-average. 

So we have to estimate the quantity X]i=i ' ^i- 
the conditions we have assumed, all terms of such a sum 
can be taken to be equal, so that we just have to estimate 
one of them. We can take • ~ — 1/| |x,i|, with |xi| ~ 
L/4 where L is the diameter of the cluster, whereas for 
estimating the modulus of the force in the direction of 
the center one can make use of ([7|). So for the velocity 
variance one gets 



with a suitable effective density poff- It is rather easily 
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where L is the hnear dimension of the cluster. In the 
case of Coma, one thus finds a value ~ 6 • lO^km^/sec^, 
which is very near to the value 5 • lO^km^/sec^ reported 
by Zwicky. 

Notice the linear dependence on L in the formula ([9]). 
In this connection one may point out that, if the external 
force were smooth, by a Taylor expansion about the ori- 
gin one would have proportional to x^, and this would 
lead to a virial (and thus also a velocity variance) propor- 
tional to rather than to L. Instead, the observations 
seem to require a proportionality to L. Apparently, this 
was first pointed out by Kazanas and Mannheim 7], in 
a paper in which some data were reported in a range of 
L covering five orders of magnitude (see 0, Fig. 2, page 
539). This property is also confirmed by a dimensional 
analysis. Indeed, with the parameters entering the prob- 
lem, the square of a velocity can be formed only as , or 
as cHqL or as (HqL)'^. But the first term is by far too 
large, the last term by far too small, while the term linear 
in L is indeed about of the correct order of magnitude. 

We finally show how the estimate (O for the effective 
density, namely, peS — 5po, is obtained. To this end 
one makes reference to the metric g^i, = 77^^ -I- h^i, with 
h^i, defined by (H]), and to the corresponding mean met- 
ric obtained by averaging with respect to the considered 
probability density for the positions of the galaxies. De- 
noting the mean by ( . ) , the mean metric is then found 
to be given by 

ds^ = ( 5^^ ) dxf'dx'' = (1 - a - 3/3) c^df ~il + a + P)df 
where dP = dx^ + dy^ + dz^ and 
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(the approximation of small velocities was used in the 
expression for /?). Introducing an effective density p^g 
such that 
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one gets 



a ~ 47rG/9offi?oVc^ , P < (2/3)a . (11) 



The consistency condition is now the requirement that 
the expansion rate calculated with the mean metric does 
actually coincide with the rate that was introduced into 
the definition of the model. This condition takes the form 
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2dt 1 - a - 3/3 
On the other hand from pTjl . using Rq — c, one gets 



(12) 
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With these expressions for a and $, the consistency con- 
dition ([T^ then becomes an algebraic one, which gives 
for peff a value that we have rounded off to 
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